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EXECUTIVE  SUMMARY 


This  report  summarizes  the  evaluation  and  refinement  of  an  environmental  stress 
index  (ESI)  for  heat  stress  assessment.  Two  independent  studies  containing  two 
different  databases  were  anaiyzed  in  order  to  evaiuate  ESI  for  different  ciimatic 
conditions  and  terrestrial  heights  beiow  and  above  sea  ievei. 

The  purpose  of  the  first  study  was  to  evaluate  the  ESI  analytically  and 
experimentally.  Meteorological  measurements  were  taken  in  different  climatic  zones 
(hot/wet,  hot/dry,  and  extremeiy  hot/dry)  at  19  different  locations  for  120  days,  and  a 
modified  ESI  was  developed  as  follows: 


ESI=0.62Ta-0.007RH+0.002SR+0.0043(Ta-RH)-0.078(0.1+SR)-^ 


Where  Ta  =  ambient  temperature  (°C),  RH  =  relative  humidity  (%),  and  SR  =  solar 
radiation  (W  m‘^).  The  correlation  coefficients  between  ESI  and  the  commonly  used  wet 
bulb  globe  temperature  (WBGT)  stress  index  were  very  high  (R^>0.898).  Therefore,  we 
conciuded  that  ESI  has  the  potential  to  be  a  practical  alternative  to  the  WBGT  based  on 
fast  response  and  accurate  climatic  microsensors  (Ta,  RH,  SR)  that  can  be  combined  in 
a  portabie  device. 

The  purpose  of  the  second  study  was  to  evaluate  the  new  and  relatively  small 
(5mm)  light  sensor  (L)  at  different  heights  beiow  and  above  sea  level  in  order  to 
measure  global  radiation  (GR)  applied  in  the  ESI  and  for  heat  stress  assessment.  Data 
were  collected  over  9  days  during  June-July  from  6  different  locations,  from  09:00 
through  17:00,  from  3  instruments:  L,  pyranometer  (P),  and  black  globe.  In  general, 
analysis  of  the  data  revealed  no  significant  differences  (P>0.05)  between  radiation 
measurements  by  P  and  L  and  between  ESI  calculated  from  L  or  P.  Therefore,  we 
concluded  that  the  L  sensor  has  the  potential  to  measure  GR  for  use  in  heat  stress 
assessment  and  in  ESI  at  different  distances  below  and  above  sea  level. 
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INTRODUCTION 


Heat  stress  evaluation  is  generally  determined  through  meteorological 
parameters  that  enable  one  to  estimate  the  influence  of  several  environmental  factors 
on  thermal  comfort  and  physiological  ability.  The  variables  included  in  heat  stress 
indices  and  their  relative  weights  have  changed  over  the  years.  Haldane  (1905) 
developed  an  index  for  heat  load  and  claimed  that  changes  in  the  wet  bulb  thermometer 
alone  were  enough  to  integrate  global  heat  load.  Other  environmental  indices  have  also 
been  suggested  which  include  measurement  of  airflow  and  thermal  radiation.  Hill  et  al. 
(10)  introduced  the  “Kata”  thermometer,  which  enabled  measurement  of  heat 
dissipation  as  a  function  of  wind  speed  and  other  parameters.  In  1932,  Vernon  (25)  was 
the  first  to  integrate  radiant  heat  into  an  environmental  stress  index  by  using  the  globe 
thermometer. 

Yaglou  and  Minard  (27)  first  introduced  the  wet  bulb  globe  temperature  (WBGT) 
during  an  extensive  period  of  development  of  heat  stress  indices  in  the  first  half  of  the 
20**^  century.  The  WBGT  index  gained  popularity  mainly  due  to  its  simplicity  and 
convenience  of  use  and  soon  was  considered  the  most  common  heat  stress  index  for 
describing  environmental  heat  stress.  This  index  is  obtained  from  three  parameters: 
black  globe  temperature  (Tg),  which  considers  the  solar  radiation;  wet  bulb  temperature 
(Tw);  and  dry  bulb  temperature  (Ta),  and  is  calculated  as  follows:  WBGT=0.7Tw  +  0.2Tg 
+  0.1  Ta.  WBGT  has  gained  immense  popularity  over  the  years  and  has  been  in  use  in 
the  field  by  the  U.S.  Army.  It  is  the  index  on  which  sports  associations  base  training 
safety  orders  as  guidance  to  prevent  heat  injury  (1,  3,  13,  14,  15).  The  World  Health 
Organization  (WHO)  has  also  adopted  it.  In  1972,  the  National  Institute  for  Occupational 
Safety  and  Health  (NIOSH)  established  the  WBGT  as  the  criterion  for  determining 
occupational  exposure  to  a  hot  environment  (18).  In  1982,  WBGT  was  approved  by  the 
ISO  as  an  international  standard  for  heat  load  assessment,  and  it  is  commonly  used  as 
a  safety  index  for  workers  in  various  occupations  (4,  5,  8,  20).  Later,  work-rest  regime 
regulations  were  made  based  on  this  index.  However,  WBGT  is  limited  in  its  evaluation 
of  heat  stress  due  to  the  inconvenience  of  measuring  Tg.  The  Tg  is  usually  measured  by 
a  thermometer  surrounded  by  a  6”  blackened  sphere,  and  purportedly  integrates  the 
global  radiation  component  of  the  thermal  load.  However,  measuring  Tg  is  cumbersome 
in  many  circumstances  for  two  main  reasons.  First,  Tg  measurement  requires  about  30 
min  for  the  instrument  to  reach  equilibrium.  Second,  the  blackened  sphere  is  often  too 
large  for  specialized  spaces  like  helicopter  cockpits  or  armored  vehicles.  Therefore, 
measuring  Tg  becomes  impractical,  especially  in  transient  situations. 

At  the  time  WBGT  was  introduced,  many  empirical  indices  were  developed  that 
originated  from  the  old  effective  temperature  (ET),  the  corrected  ET  (GET),  the  modified 
ET  (MET),  and  the  equivalent  ET  (ETR)  (2,  11,  26).  All  of  these  indices  were  based  on 
Ta,  Tw,  andTg,  but  none  was  ever  associated  with  any  physiological  variable.  The  WBGT 
has  been  only  partly  evaluated  by  analyzing  the  association  between  heat  stroke  cases 
in  soldiers  and  the  heat  load  assessed  by  the  WBGT  (27).  In  fact,  there  is  no  known 
laboratory-  controlled  study  where  the  WBGT  was  evaluated.  Recently,  Moran  et  al. 
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(16)  introduced  a  new  environmental  stress  index  (ESI)  based  on  measurements  of  Ta, 
relative  humidity  (RH),  and  solar  radiation  (SR).  ESI  was  highly  correlated  with  the 
WBGT  index  (16).  The  ESI  as  a  stress  index,  for  the  first  time,  incorporates  direct 
measurements  of  SR  and  RH.  Additionally,  all  the  three  meteorological  variables  that 
make  up  the  ESI  are  characterized  by  fast  reading  sensors  easily  obtained 
commercially.  However,  ESI  needs  to  be  further  evaluated  and  validated  under  different 
climatic  and  topographical  conditions. 

The  flow  of  global  radiation  falling  to  the  earth’s  surface  is  dependent  on  four 
factors.  The  first  factor  is  astronomical,  influenced  by  the  distance  of  the  sun  from  the 
measured  site  and  the  length  of  daylight.  The  distance  of  the  sun  from  a  site  is 
dependent  on  the  topographical  height  of  the  surface  and  the  hour  of  the  day.  With 
higher  topographical  height,  thermal  radiation  becomes  more  intense.  Thermal  radiation 
is  also  more  intense  as  daylight  becomes  extended.  The  second  factor  is  the 
transparency  of  the  atmosphere.  This  is  a  measurement,  which  changes  daily  and 
monthly.  Radiation  is  more  intense  as  atmosphere  becomes  clearer.  The  third  important 
factor  is  the  amount  of  clouds  in  the  sky.  Cloudy  skies  greatly  influence  the  measured 
global  radiation.  With  cloudy  skies,  thermal  radiation  becomes  weaker.  Cloudiness,  of 
course,  changes  as  the  seasons  change  during  the  year,  but  it  can  also  change  daily. 
The  fourth  and  final  factor  is  the  surface  coverage  (Albedo).  Different  types  of  surfaces 
are  influenced  in  variable  ways  depending  on  the  reflected  solar  radiation.  The  more 
light  is  reflected  from  an  exterior  surface,  the  greater  the  measured  radiation.  A  snow- 
covered  surface  is  a  good  example  of  a  surface  greatly  influenced  by  the  reflected 
radiation.  Here  the  measured  radiation  is  greater  than  a  blackened  surface  that  absorbs 
the  sun’s  radiation.  Therefore,  any  device  required  to  measure  solar  radiation  should 
also  be  validated  by  all  these  factors. 

Over  the  last  few  years,  major  progress  has  been  made  in  the  field  of  portable 
heat  load  instrument  technology.  New  instruments  are  able  to  measure  and  present  a 
variety  of  meteorological  parameters  (e.g.,  Ta,  RH,  and  wind  speed)  and  display  the 
calculated  heat  load  on  a  screen.  Moreover,  no  special  training  is  needed  to  operate 
these  instruments.  It  is  worth  mentioning  that  the  microsensors  in  these  instruments 
have  a  fast  response  time,  and  displaying  the  results  is  a  matter  of  waiting  just  a  few 
seconds.  The  simplicity  of  operation  makes  these  instruments  previously  inaccessible  to 
laymen. 

There  were  two  purposes  for  the  present  study.  The  first  was  to  evaluate  and 
validate,  using  a  large  database,  the  thermal  load  index  (ESI)  constructed  from 
microsensors  for  Ta,  RH  and  solar  radiation  (SR),  and  to  assess  and  determine  whether 
this  newly  developed  index  can  serve  as  a  reliable  and  valid  alternative  to  WBGT  for 
measuring  environmental  stress.  The  second  purpose  was  to  evaluate  a  new  and 
relatively  small  (5  mm)  light  sensor  (L)  that  is  easily  deployable  at  different  heights 
above  sea  level  in  order  to  measure  global  radiation  in  the  ESI. 
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MILITARY  RELEVANCE 


It  is  important  for  military  commanders  to  have  knowledge  in  real-time  about 
environmental  stress  occurring  during  different  scenarios  (e.g.,  during  training,  combat 
and  deployments  to  hot/cold  environments).  Current  techniques  of  evaluating 
environmental  stress  are  not  facile.  They  are  also  cumbersome,  labor  intensive,  time 
consuming,  and  do  not  use  modern  technology  or  advanced  sensor  systems  for 
weather  measurements.  Currect  systems  with  reliable  microsensors  also  have  not  been 
fully  integrated  with  microprocessor  units  into  a  portable,  ultra  lightweight,  friendly  to  use 
device. 

The  existing  U.S.  military  heat  stress  monitoring  systems  are  based  largely  on 
the  WBGT  (14),  whereas  the  IDF  system  is  based  on  the  Dl  (22).  The  inherent 
limitations  of  the  WBGT  and  the  Dl  have  been  reported  (7,12,15)  in  terms  of 
applicability  across  a  broad  range  of  potential  military  scenarios  and  environments. 
These  limitations  are  attributed,  in  part,  to  early  constraints  on  sensor  and 
computational  complexity. 

Mathematical  models  of  human  heat  strain  allow  full  consideration  of  the  complex 
interactions  of  environment,  clothing,  acclimatization  status,  and  metabolic  heat 
production  that  ultimately  determine  a  soldier’s  performance  limits  in  a  given  scenario. 


METHODS 


STUDY  I 

The  database  was  obtained  from  a  study  collected  at  various  climatic  sites  in 
Israel  and  contained  126,558  measurements  for  each  of  the  variables  measured. 
Weather  measurements  were  collected  every  10  min,  during  24  hours  a  day,  for  120 
days,  at  19  different  locations  around  Israel,  including  sites  near  the  Mediterranean  sea 
characterized  as  a  high  humidity  zone,  sites  in  the  desert  characterized  as  a  hot/dry 
zone,  and  sites  near  the  Red  Sea  characterized  as  an  extremely  hot/dry  zone. 
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Measurements 


The  official  Israeli  Meteorological  Service  collected  weather  measurements  for 
the  nineteen  locations.  Ta  and  Tw  were  measured  with  Campbell  thermometers  (model 
HMP45C),  and  relative  humidity  (RH)  was  measured  with  a  Rotronic  instrument  (model, 
MP  100A).  These  three  instruments  were  placed  under  a  shelter  (Stevenson  screen). 
Under  open  sky,  Tg  was  measured  using  the  Vernon  black  globe  thermometer,  solar 
radiation  (SR)  was  measured  using  the  EPLAB  radiometer  (sensitivity  of  285-2800  nm), 
and  the  infra-red  (IR)  light  sensor  measured  by  Centro  Vision  sensor  (model  CD-1705) 
with  peak  sensitivity  of  800-920  nm. 

Calculations 


Heat  stress  indices  were  calculated  as  follows:  WBGT  was  calculated  according 
to  Yaglou  and  Minard’s  standard  formula;  WBGT  =  0.7Tw  -i-0.2Tg  +0.1  Tg  (22).  The  newly 
modified  environmental  stress  index  (ESI)  was  calculated  as  follows: 

ESI=0.62Ta-0.007RH+0.002SR+0.0043(TaRH)-0.078(0.1+SR)‘^ 


Statistical  Analysis 

For  the  evaluation  of  the  modified  ESI,  we  validated  a  series  of  models  for  WBGT 
as  dependent  variables,  and  Ta,  RH,  SR,  their  interactions  and  different  transformations 
(inverse,  and  quadratic)  as  independent  variables.  Pearson’s  analysis  was  used  for 
testing  bivariate  correlations  between  the  independent  and  the  dependent  variables.  All 
models  were  linear  models  and  used  the  least  squares  algorithm.  Optimization  of  the 
constants  was  executed  by  the  DUD  (does  not  use  derivative)  method  (19).  For  all 
models,  we  computed  the  coefficient  of  determination  (R-Square  of  the  model)  and 
plotted  a  series  of  residual  plots  versus  predicted  values  for  all  data  and  for  every 
meteorological  station  separately.  All  statistical  contrasts  were  accepted  at  the  P<0.05 
level  of  significance.  Data  are  presented  in  this  study  as  meanslSE.  For  all 
computations  and  statistical  analysis,  we  used  SAS  8.0  Software,  Procedures  CORR 
and  GLM. 


STUDY  II 

This  database  was  obtained  from  a  study  conducted  in  Israel,  during  the  summer 
season  (June-July)  for  9  days  at  6  different  meteorological  stations  at  different  terrestrial 
heights,  from  the  lowest  place  in  the  world  at  -400m,  located  at  the  Dead  Sea,  and  up  to 
1600m,  located  at  a  ski  resort  in  Mt.  Hermon.  Global  or  solar  radiation  measurements 
were  collected  for  15-min  daily  between  09:00  and  17:00.  These  data  were  collected 
from  the  following  three  instruments  placed  under  the  open  sky  at  1  m  height:  black 
globe  (Tg)  using  the  Vernon-black  globe  thermometer,  EPLAB  pyranometer  (P),  and 
infra-red  light  (IRL)  sensor. 
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Statistical  Analysis 


The  new  model  to  predict  P  values  from  L  data  was  validated  as  a  linear  model 
that  was  fitted  by  the  least  square  method  using  SAS  8.0  software.  Correlation 
coefficients  between  ESI  calculated  from  P  (ESIp)  and  IRL  sensor  (ESIir)  were 
computed  using  Pearson  correlation  analysis.  All  statistical  contrasts  were  accepted  at 
the  P<0.01  level  of  significance. 


RESULTS 


STUDY  I 

These  data  were  collected  every  15  min  over  24  hr  for  120  days.  Therefore,  a 
wide  range  of  weather  measurements,  over  1 26,000  for  each  variable,  was  covered  as 
depicted  in  Table  1  (mean  ±  SD),  and  for  each  of  the  19  stations  in  Appendix  A. 


Table  1.  Mean  (±  SD)  and  range  of  environmental  measurements  of  the  ESI  validation 
versus  the  WBGT.  Data  were  collected  from  1 9  different  locations  every  1 5  min  over  24 
hours  for  120  days. 


Ta 

Tw 

RH 

Tg 

SR 

WBGT 

(°C) 

(°C) 

(%) 

(°C) . . 

(Wm-2) 

(°C) 

Mean±SD 

23.0716.35 

16.8813.66 

55.52124.00 

26.40110.81 

2891353 

19.0514.71 

Range 

0.80-44.30 

0.60-29.70 

1.75-100 

0-59.42 

0-1337 

0.7-32.58 

In  order  to  better  correlate  ESI  with  WBGT  and  to  improve  the  distribution  around 
the  line  of  identity,  we  modified  the  ESI  constants.  This  newly  modified  ESI  is  based  on 
the  same  three  variables  (Ta,  RH,  and  SR),  interaction  (Ta-RH)  and  transformation 
[(SR+0.1)'^]  as  the  original  ESI  (16).  The  ESI  modification  is  a  result  of  statistical 
analysis,  which  included  three  steps.  The  first  step  applied  the  original  ESI  on  the  19 
different  databases  and  analyzed  them  for  correlation  between  WBGT  and  ESI,  and  the 
residual  scattergram  from  the  line  of  identity.  In  the  second  step,  we  validated  the 
WBGT  as  a  dependent  variable,  and  each  of  the  variables  and  components  in  the 
original  ESI  as  an  independent  variable.  Since  the  results  from  the  second  step 
revealed  a  very  high  correlation  between  the  dependent  and  the  independent  variables, 
we  progressed  to  the  third  step.  At  the  second  stage,  the  ESI  concept  and  construction 
was  evaluated.  For  the  third  step,  we  optimized  the  different  constants  of  the  ESI 
variables,  interaction  and  transformation  as  follows: 

For  Ta  measured  In  °C: 
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ESI=0.62Ta-0.007RH+0.002SR+0.0043(TaRH)-0.078(0.1+SR)'^ 

For  Ta  measured  in  °F: 

ESI=0.73Ta-0.007RH+0.005SR+0.0026(TaRH)-0.11 5(0.1 +SR)‘^ 

The  modified  ESI  was  applied  separately  to  each  of  the  19  databases,  and  overall,  a 
highly  significant  correlation  coefficient  (R^>0.898,  P<0.001)  was  obtained  with  residuals 
distributed  symmetrically  around  the  zero  line  in  most  of  the  locations  as  depicted  in 
Figures1-19. 
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ESI  (units)  Residuals 


Figure  1 .  Comparison  of  the  modified  ESI  with  the  WBGT  index  showing 
correlation  (bottom)  and  residuals  scattergram  (top).  Database  for 
this  figure  was  collected  from  station  1 
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Figure  2.  Comparison  of  the  modified  ESI  with  the  WBGT  index  showing 
correlation  (bottom)  and  residuals  scattergram  (top).  Database  for 
this  figure  was  collected  from  station  2 


Figure  3.  Comparison  of  the  modified  ESI  with  the  WBGT  index  showing 
correlation  (bottom)  and  residuals  scattergram  (top).  Database  for 
this  figure  was  collected  from  station  3 
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Figure  4.  Comparison  of  the  modified  ESI  with  the  WBGT  index  showing 
correlation  (bottom)  and  residuals  scattergram  (top).  Database  for 
this  figure  was  collected  from  station  4 
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ESI  (units)  Residuals 


Figure  5.  Comparison  of  the  modified  ESI  with  the  WBGT  index  showing 
correlation  (bottom)  and  residuals  scattergram  (top).  Database  for 
this  figure  was  collected  from  station  5 
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Figure  6.  Comparison  of  the  modified  ESI  with  the  WBGT  index  showing 
correlation  (bottom)  and  residuals  scattergram  (top).  Database  for 
this  figure  was  collected  from  station  6 
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Figure  7.  Comparison  of  the  modified  ESI  with  the  WBGT  index  showing 
correlation  (bottom)  and  residuals  scattergram  (top).  Database  for 
this  figure  was  collected  from  station  7 
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Figure  8.  Comparison  of  the  modified  ESI  with  the  WBGT  index  showing 
correlation  (bottom)  and  residuals  scattergram  (top).  Database  for 
this  figure  was  collected  from  station  8 
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Figure  9.  Comparison  of  the  modified  ESI  with  the  WBGT  index  showing 
correlation  (bottom)  and  residuals  scattergram  (top).  Database  for 
this  figure  was  collected  from  station  9 
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ESI  (units)  Residuals 


Figure  10.  Comparison  of  the  modified  ESI  with  the  WBGT  index  showing 
correlation  (bottom)  and  residuals  scattergram  (top).  Database  for 
this  figure  was  collected  from  station  10 


WBGT  (“C) 
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ESI  (units)  Residuals 


Figure  1 1 .  Comparison  of  the  modified  ESI  with  the  WBGT  index  showing 
correlation  (bottom)  and  residuals  scattergram  (top).  Database  for 
this  figure  was  collected  from  station  1 1 


WBGT  (“C) 
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ESI  (units)  Residuals 


Figure  12.  Comparison  of  the  modified  ESI  with  the  WBGT  index  showing 
correlation  (bottom)  and  residuals  scattergram  (top).  Database  for 
this  figure  was  collected  from  station  12 
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Figure  13.  Comparison  of  the  modified  ESI  with  the  WBGT  index  showing 
correlation  (bottom)  and  residuals  scattergram  (top).  Database  for 
this  figure  was  collected  from  station  1 3 


ESI  (units)  Residuals 


Figure  14.  Comparison  of  the  modified  ESI  with  the  WBGT  index  showing 
correlation  (bottom)  and  residuals  scattergram  (top).  Database  for 
this  figure  was  collected  from  station  14 
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Figure  15.  Comparison  of  the  modified  ESI  with  the  WBGT  index  showing 
correlation  (bottom)  and  residuals  scattergram  (top).  Database  for 
this  figure  was  collected  from  station  15 


Figure  16.  Comparison  of  the  modified  ESI  with  the  WBGT  index  showing 
correlation  (bottom)  and  residuals  scattergram  (top).  Database  for 
this  figure  was  collected  from  station  16 


ESI  (units)  Residuals 


Figure  17.  Comparison  of  the  modified  ESI  with  the  WBGT  index  showing 
correlation  (bottom)  and  residuals  scattergram  (top).  Database  for 
this  figure  was  collected  from  station  1 7 


Figure  18.  Comparison  of  the  modified  ESI  with  the  WBGT  index  showing 
correlation  (bottom)  and  residuals  scattergram  (top).  Database  for 
this  figure  was  collected  from  station  18 


ESI  (units)  Residuals 


Figure  19.  Comparison  of  the  modified  ESI  with  the  WBGT  index  showing 
correlation  (bottom)  and  residuals  scattergram  (top).  Database  for  this  figure  was 
collected  from  station  1 9 


STUDY  II 


The  dataset  collected  in  this  study  included  global  radiation  measured  by 
pyranometer  (P),  infra-red  (IR)  light  sensor,  and  black  globe  thermometer  (Tg),  wet  bulb 
temperature  (Tw),  and  ambient  temperature  (Ta).  The  measurements  were  carried  out 
over  a  wide  range  of  weather  measurements  as  depicted  in  Table  2  for  meanlSD 
values  and  in  Appendix  B  for  each  of  the  6  locations. 

There  were  no  significant  differences  (NS,  P>0.05)  between  the  measurements 
of  the  three  IR  sensors  as  depicted  in  Table  2  and  for  each  of  the  6  locations  in 
Appendix  B.  However,  higher  values  for  the  pyranometer  in  comparison  to  the  3  IR 
sensors  measurements  were  found  mainly  between  11:00  and  14:00  at  the  6  different 
locations  (see  Figures  20-25).  Analysis  of  the  Tg  measurements  showed  similar  patterns 
as  in  the  radiation  sensors  (P  and  L)  measurements. 


Table  2.  Mean  (+SD)  and  range  of  global  radiation  (GR)  measured  using  2 
pyranometers  [by  Heller  Institute]  (P)  and  by  Israeli  Meteorological  Service  (Pims)].  black 
globe  thermometer  (Tg),  and  3  infra-red  (IR)  light  sensor.  Database  was  collected  from  6 
different  meteorological  stations  at  various  distances  below  (-200  m  and  -400  m)  and 
above  sea  level  (n=238). 


P 

P|MS 

Tg 

IRi 

IR2 

IR3 

Tw 

(Wm-2) 

(Wm‘^) 

(°C) 

(Wm'^) 

(Wm'^) 

(Wm'2) 

(°C) 

(°C) 

Mean 

±SD 

754± 

189 

81 8± 
254 

42.35± 

4.76 

651± 

172 

639± 

158 

665± 

157 

28.69± 

4.39 

21 .52± 
2.14 

Range 

205-1028 

40-1058 

30-52 

151-875 

160-876 

179-887 

19.1-38.6 

16.5-26.0 
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Figure  20.  Global  radiation  (GR)  measured  by  2  pyranometers  (P,  Pims),  black 
globe  thermometer  (Tg),  and  3  infra-red  (IR)  light  sensors  at  -400  m  below  sea  level 
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Figure  21 :  Global  radiation  (GR)  measured  by  pyranometer  (P),  3  infra-red  (IR) 

light  sensors  and  black  globe  thermometer  (Tg)  at  -200  m  below  sea  level. 
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Figure  23:  Global  radiation  (GR)  measured  by  pyranometer  (P|ms)>  irrfra-red  (IR) 

light  sensor  and  black  globe  thermometer  (Tg)  at  400  m  from  sea  level. 
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The  ESI  was  calculated  separately  using  the  three  IR  sensors  (ESIim,  ESIir2, 
ESI1R3)  and  the  pyranometer  (ESIp)  measurements,  and  applied  separately  to  each  of 
the  6  locations  differing  in  their  distances  below  (-200  m  and  -400  m)  or  above  sea  level 
(Figures  26-31).  No  significant  differences  were  found  between  the  3  different  ESIs 
calculated  with  IR  sensor  or  the  pyranometer  and  the  WBGT,  as  depicted  in  Table  3, 


Figure  25:  Global  radiation  (GR)  measured  by  pyranometer  (P),  2  infra-red  (IR)  light 
sensors  and  black  globe  thermometer  (Tg)  at  1600  m  from  sea  level. 
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and  for  each  of  the  locations  in  Appendix  C.  However,  when  the  ESI  was  compared  to 
the  WBGT  index,  no  significant  differences  were  found  in  4  locations  (-400  m  below  and 
30  m,  400  m  and  900  m  above  sea  level).  At  -200m  below  sea  level,  WBGT  was 
significantly  lower  (P<0.05)  than  ESI,  and  at  1600m  above  sea  level,  WBGT  values 
were  significantly  higher  (P<0.05)  than  the  ESI  values. 
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Table  3.  Mean  (±SD)  and  range  of  the  modified  ESI  and  the  WBGT  index  calculated 
from  global  radiation  measured  by  pyranometer  (ESIp)  and  3  different  IRL 
sensors  (ESIir).  Database  was  collected  from  6  different  meteorological  stations  at 
different  heights  below  and  above  sea  level  (n=238). 


ESIp 

ESIim 

ESIir2 

ESIir3 

WBGT(°C) 

Mean±SD 

26.03+3.11 

26.31  ±3.00 

26.50±2.88 

27.21  ±2.76 

26.39±2.63 

Range 

18.68-32.63 

19.28-32.21 

19.79-32.29 

20.61-32.33 

18.37-30.76 
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Figure  26:  Comparison  of  the  modified  ESI  with  the  WBGT  index.  ESI  was  calculated  from  global 
radiation  (GR)  measured  by  pyranometer  (ESip)  and  3  different  IR  light  sensors 
(ESI|p)  at  -  400  m  below  sea  level. 


Figure  28:  Comparison  of  the  modified  ESI  with  the  WBGT  index.  ESI  was  calculated  from  GR  measured 
by  pyranometer  (ESlp)  and  3  different  iR  light  sensors  (ESI|p)  at  30  m  from  sea  level. 
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Figure  29:  Comparison  of  the  modified  ESI  with  the  WBGT  index.  ESI  was  calculated  from  GR  measured 
by  pyranometer  (ESlp)  and  IR  light  sensors  (ESI,p)  at  400  m  from  sea  level. 


Figure  30:  Comparison  of  the  modified  ESI  with  the  WBGT  index.  ESI  was  calculated  from  GR  measured 
by  pyranometer  (ESIp)  and  2  different  IR  light  sensors  (ESIjp)  at  900  m  from  sea  level. 


DISCUSSION 


STUDY  I 

The  modified  ESI  introduced  in  this  study  to  evaluate  environmental  heat  stress 
reliably  matches  values  calculated  by  using  the  conventional  WBGT.  This  modified  ESI 
is  based  on  the  same  parameters  as  ESI  (Ta,  RH  and  SR),  which  integrates  the  thermal 
load  derived  by  the  specific  climatic  conditions.  However,  separate  analysis  of  the  1 9 
different  databases  obtained  from  the  1 9  different  locations  revealed  changes  in  the  ESI 
constants  of  the  variables,  interaction,  and  transformation.  Arithmetically,  the  first 
developed  ESI  was  well  correlated  to  WBGT  (16).  However,  validation  of  the  modified 
ESI  for  these  large  databases  further  improved  the  correlation  with  the  WBGT 
((R^>0.898).  The  strength  of  any  prediction  index  is  its  ability  to  predict  with  high 
correlation  the  measured  or  calculated  universal  index  under  widely  varied  climatic 
conditions. 

ESI  differs  from  other  indices  that  have  been  suggested  in  the  past  in  two  critical 
ways.  First,  this  stress  index  for  the  first  time  uses  direct  measurements  of  SR  and  RH. 
These  direct  measurements  of  RH  and  SR,  when  used  in  ESI,  are  not  as  cumbersome 
as  measuring  Tg  and  Tw  for  calculating  the  WBGT.  Second,  the  three  meteorological 
variables  used  in  ESI  are  characterized  by  fast-reading  responses  that  take  only  a  few 
seconds  to  reach  equilibrium.  For  an  index  to  be  valid  and  practical,  it  should  allow  the 
comparison  and  evaluation  of  a  combination  of  different  meteorological  parameters  as 
far  as  their  influence  on  the  individual  is  concerned.  It  also  helps  to  find  different 
combinations  of  these  parameters  that  cause  equal  subjective  heat  sensations. 
Moreover,  the  index  must  enable  one  to  assess  the  different  weights  of  each  of  the 
meteorological  parameters  on  the  individual  (6). 

In  conclusion,  this  study  separately  confirms  the  results  of  a  previous  study  (16), 
which  initially  introduced  ESI  and  found  a  high  correlation  with  WBGT.  However,  in  the 
present  study,  ESI  was  evaluated  using  a  larger  database  with  revised  constants 
included  in  its  determination  and  a  higher  correlation  with  WBGT.  Therefore,  the  ESI 
has  the  potential  to  serve  as  a  substitute  to  the  WBGT  index  and  to  be  used  for  safety 
limits  during  training  and  military  activity.  It  can  also  serve  as  a  part  of  the  guidelines  for 
work  rest  cycles  and  fluid  replacement.  To  implement  the  current  guidelines  and 
limitations  for  exercise  in  a  hot  climate,  there  is  a  great  need  for  the  development  of  an 
accurate,  portable  heat  stress  measurement  device.  This  study  suggests  that  the  ESI  is 
an  ideal  index  to  be  used  in  such  a  device.  However,  further  evaluation  of  ESI  is 
required  for  validation  with  specific  physiological  variables  (e.g.,  sweating  rate,  heart 
rate,  and  internal  body  temperatures). 
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STUDY  II 

In  a  previous  study  (16),  a  newly  developed  algorithm  from  light  sensor 
appeared  valid  for  use  in  determining  predictive  global  radiation  (PGR),  since  the 
correlation  coefficient  between  PGR  and  measured  global  radiation  by  pyranometer  was 
very  high.  In  this  study,  we  additionally  evaluated  PGR  at  different  heights  below  and 
above  sea  level,  and  the  database  for  this  study  was  collected  at  6  separate  locations. 
The  range  of  the  measurements  for  each  day  was  wide  as  seen  in  Table  2.  Analyses  of 
the  global  radiation  measurements  at  different  heights  above  and  below  sea  level  is 
limited  in  this  study  because  a  real  comparison  can  be  done  only  from  a  large  database 
containing  annual  GR  measurements  collected  at  the  same  time  at  the  6  locations.  In 
an  analysis  made  at  different  heights  from  sea  level,  which  were  made  by  the  IMS  over 
26  years,  it  was  found  that  at  -208  m  below  sea  level,  GR  values  were  lower  by  100 
W  m’^  from  measurements  at  815  m  above  sea  level.  These  differences  are  consistent, 
but  vary  according  to  the  month  of  the  year.  However,  mean  GR  values  were  higher  at 
-395m  than  at  -208m.  The  latter  might  be  explained  by  other  factors  that  influence  the 
GR  apart  from  the  topographical  height  from  sea  level  (e.g.,  cloudiness,  atmospheric 
transparency,  and  reflective  radiation). 

In  analyzing  the  data  from  the  3  IR  light  sensors  and  the  pyranometer,  we  can 
draw  two  conclusions.  First,  there  were  no  significant  differences  between  these  3 
sensors,  which  strengthens  the  reliability  and  validity  of  the  IR  light  sensor  to  measure 
GR.  The  second  conclusion  relates  to  a  comparison  between  the  P  measurements  and 
the  IR  light  sensors.  Higher  values  for  P  were  measured  in  all  the  6  locations.  Although 
differences  were  not  significant  for  all  the  locations,  for  better  accuracy,  refinement  of 
the  model  that  predicts  global  radiation  (PGR)  from  the  IR  light  sensor  should  be 
considered. 

In  this  study,  ESI  was  applied  by  using  3  different  IR  light  sensors  and 
pyranometer  measurements  for  global  radiation.  However,  no  statistically  significant 
differences  were  found  between  the  ESI  calculated  with  either  sensor  (Figures  26-31). 
Overall,  when  WBGT  was  applied  to  the  same  6  databases,  no  statistically  significant 
differences  were  found  in  4  databases  between  ESI  and  WBGT.  However,  at  -200  m 
below  sea  level,  WBGT  values  were  significantly  lower  (P  <0.05)  than  ESI,  and  at  1600 
m  above  sea  level,  WBGT  values  were  higher  than  ESI. 

The  present  study  derived  a  simple  fast  response  method  to  measure  and 
evaluate  global  radiation.  This  method  is  easier  to  measure,  interpret,  and  use  as  a 
potential  substitute  for  Tg,  which  is  incorporated  in  the  WBGT. 
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CONCLUSIONS 


Although  the  WBGT  is  the  most  common  heat  stress  index  and  has  been 
adopted  by  many  organizations,  we  introduced  and  evaluated  in  this  study  a  new 
environmental  stress  index  (ESI)  that  may  be  a  reasonable  substitute  for  the 
conventional  index.  The  ESI  was  fabricated  using  new  micro-sensors  that  measure 
climatic  conditions  with  programmable  microchips.  The  present  study  confirms  that  the 
simple  ESI  is  valid  to  evaluate  wide  climatic  conditions.  This  simple  index  should  be 
easier  to  measure  since  it  is  based  on  more  commonly  used  sensors  (e.g.,  RH,  GR) 
than  the  WBGT  index.  The  ESI  is  capable  of  overcoming  the  technical  difficulties  of 
measuring  WBGT  (e.g.,  Tg,  Tw),  while  providing  an  easy  to  operate  index  based  on  fast 
response  micro-sensors  to  assess  heat  stress.  This  development  is  very  encouraging 
since  there  already  exist  micro-sensors  for  measuring  Ta  and  relative  humidity  (RH)  that 
are  available  off-the-shelf.  We  suggest  that  ESI  has  the  potential  to  be  widely  accepted, 
used  universally,  and  easily  implemented  in  a  relatively  small  portable  device.  However, 
further  investigation  is  necessary  in  order  to  refine  the  PGR  distinct  from  the  IR-light 
sensor. 
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APPENDIX  A.  Mean  (±  SD)  and  range  of  environmental  measurements  of  the  ESI 
validation  versus  WBGT  at  the  1 9  different  locations 


1  o 

Tw 

(°C) 

RH 

(%) 

Tg 

(°C) 

SR 

(Wm’2) 

WBGT 

(°C) 

Station  no.  1 

19.29±5.75 

3.3612.51 

53125 

22.12110.23 

2771352 

15.7113.96 

8.10-34.80 

6.30-21.80 

6-98 

5.00-50.00 

0-1172 

6.87-26.21 

Station  no.  2 

20.33±4.84 

15.9512.48 

65119 

25.4819.  08 

2401308 

18.3013.70 

8.90-39.20 

7.90-22.60 

5-67 

8.07-50.01 

0-993 

8.15-27.94 

Station  no.  3 

19.35±5.29 

15.3212.20 

68123 

22.2419.97 

2311296 

17.1113.75 

9.50-40.60 

9.10-22.00 

9-100 

6.40-49.20 

0-972 

8.69-27.88 

Station  no.  4 

26.0014.73 

16.2812.00 

35116 

28.2918.71 

2801354 

19.6513.14 

15.40-39.00 

10.40-21.50 

5-92 

8.70-53.20 

0-1107 

10.88-28.00 

Station  no.  5 

19.9515.97 

15.3713.30 

64124 

19.76110.97 

2591325 

16.7114.67 

5.30-37.80 

4.90-29.70 

10-100 

5.10-45.70 

0-1143 

4.28-31.67 

Station  no.  6 

25.2015.21 

16.4112.10 

39117 

30.09110.58 

2471314 

20.0213.63 

14.20-38.60 

8.40-22.70 

5-95 

11.86-59.42 

0-989 

10.03-30.26 

Station  no.  7 

17.7615.22 

11.7912.28 

50126 

21.8319.94 

2821363 

14.4013.59 

9.30-32.40 

5.40-19.70 

5-97 

5.86-50.70 

0-1337 

6.31-25.88 

Station  no.  8 

15.3815.75 

12.2313.49 

72123 

17.6019.51 

2701343 

13.6214.73 

0.80-33.20 

0.60-22.60 

26-100 

0-47.10 

0-1152 

0.70-28.28 

Station  no.  9 

27.5913.96 

18.3911.81 

39111 

33.6418.27 

2661338 

22.3612.80 

18.40-38.80 

12.30-25.40 

10-99 

12.56-59.33 

0-1089 

14.66-29.84 

Station  no.  1 0 

20.8315.61 

15.8012.89 

61120 

23.93110.56 

2621331 

17.9314.37 

9.20-37.00 

8.00-24.00 

14-94 

5.92-52.80 

0-1131 

7.81-30.30 

Station  no.  1 1 

21.0514.12 

16.4811.95 

64120 

20.69110.18 

2771338 

17.7813.34 

13.90-36.40 

10.30-23.00 

14-95 

10.90-54.20 

0-1055 

9.27-29.73 

Station  no.  12 

24.3715.89 

17.5311.77 

54126 

28.18111.11 

3341385 

20.3413.65 

13.70-37.70 

12.00-25.10 

10-97 

11.20-56.10 

0-1063 

12.35-29.54 

Station  no.  13 

24.9913.61 

20.7811.96 

68115 

31.2318.65 

3121367 

23.2913.14 

15.10-33.30 

13.90-24.80 

27-97 

14.50-46.60 

0-1098 

14.22-29.19 

Station  no.  14 

24.3514.37 

20.5011.91 

72120 

28.58110.56 

3051355 

22.5113.53 

14.90-34.60 

13.70-25.70 

25-100 

12.50-50.40 

0-1073 

13.58-30.45 

Station  no.  1 5 

32.6414.75 

19.5111.46 

27113 

35.9318.26 

3281379 

24.1112.72 

20.80-44.30 

15.50-25.00 

2-68 

19.00-55.00 

0-1038 

16.94-30.69 

Station  no.  1 6 

25.8815.09 

20.3612.21 

62120 

25.53110.69 

3241367 

21.9513.77 

13.90-40.80 

10.80-25.60 

10-100 

5.20-52.70 

0-1022 

10.07-30.54 
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Station  no.  17 

Station  no.  18 

Station  no.  19 


23.4814.49 

15.7512.15 

44122 

28.61110.03 

3491403 

13.90-36.30 

9.10-22.30 

2-96 

12.00-54.70 

0-1147 

26.3814.65 

21.0312.40 

63119 

29.9919.84 

3201367 

15.50-41.60 

12.60-27.80 

17-94 

12.10-56.60 

0-1038 

24.0012.67 

19.7611.99 

67±17 

24.9019.74 

3181366 

18.10-33.40 

12.20-25.10 

16-95 

9.00-46.50 

0-1034 

1 9.0913.38 
11.84-29.03 

23.3613.69 

13.16-32.58 

21.2112.94 

13.81-29.59 
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APPENDIX  B.  Mean  (±SD)  and  range  of  global  radiation  (GR)  measured  using  2 
pyranometers  [by  Heller  Institute  (P)  and  by  Israeli  Meteorological  Service  (Pims)],  black 
globe  thermometer  (Tg),  and  3  infra-red  (IR)  light  sensors.  Database  was  collected  from 
6  different  meteorological  stations  at  different  distances  (Dist.)  above  and  below  sea 

level. 


Dist. 

from  SI 

(m) 

P 

(Wm'^) 

P|MS 

(Wm-2) 

Tg 

(°C) 

IRi 

(Wm'^) 

IR2 

(Wm-2) 

IR3 

(Wm’^) 

Ta 

(°C) 

Tw 

(°C) 

-400 

8171180 

917177 

49.2011.52 

7211127 

6721152 

7351152 

35.7512.05 

22.8210.85 

387- 

722-992 

45-52 

410-849 

325-828 

370-887 

32.20- 

22.00- 

1028 

38.60 

24.50 

-200 

7481214 

40.7114.45 

5811204 

6371161 

6741150 

31 .9912.99 

24.1710.96 

205-967 

30-48 

151-821 

160-814 

179-844 

26.10- 

22.50- 

36.40 

26.00 

30 

7071183 

8101187 

40.1112.01 

6451165 

6281181 

6831179 

27.4711.12 

22.3110.72 

290- 

368-990 

35-45 

296-819 

192-840 

220-871 

24.20- 

21 .00- 

1028 

28.90 

23.00 

400 

7091348 

6881155 

6941139 

26.1114.16 

20.6811.42 

40-1053 

237-875 

292-876 

19.10- 

18.20- 

31.10 

22.50 

900 

8061159 

9131153 

38.0312.20 

24.9411 .56 

19.0311.31 

471- 

560- 

32-40.50 

21.90- 

16.50- 

1003 

4058 

27.90 

21.00 

1600 

7081177 

45.0712.63 

5721133 

5791105 

27.5010.86 

20.1111.84 

411-931 

39-48.50 

301-724 

346-704 

25.50- 

17.00- 

28.90 

22.50 
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APPENDIX  C.  Mean  {±SD)  and  range  of  the  modified  ESI  and  the  WBGT  index 
calculated  from  global  radiation  measured  by  pyranometer  (ESIp)  and  3  different 
infra-red  light  sensors  (ESIir).  Database  was  collected  from  6  different  meteorological 
stations  at  different  distances  (Dist.)  above  and  below  sea  level. 


Dist.  from  SL 
(m) 

ESIp 

ESIim 

ESI|R2 

CO 

DC 

CO 

LU 

WBGT 

(°C) 

-400 

29.59±0.77 

29.19-31.10 

29.4010.86 

27.89-30.93 

29.3010.77 

27.90-30.70 

29.4210.77 

28.02-30.83 

29.3910.86 

27.67-30.76 

-200 

29.60±1.69 

25.01-32.63 

29.2611.67 

24.90-32.21 

29.3811 .61 
24.92-32.29 

29.4511.57 

24.96-32.33 

28.2611.32 

24.71-30.25 

30 

26.21±1.12 

23.32-28.18 

26.0911 .08 
23.31-27.74 

26.0511.19 

22.99-27.76 

26.1611.19 

23.07-27.81 

26.3810.85 

24.23-27.52 

400 

23.7812.71 

18.37-26.60 

23.7112.68 

18.68-26.57 

23.7812.71 

18.37-26.60 

900 

24.0611.62 

19.90-26.47 

23.8211.58 

19.75-26.40 

23.8311.63 

19.79-26.24 

23.4211.43 

20.14-25.36 

1600 

23.9611.55 

20.63-25.85 

23.6911.49 

20.51-25.59 

23.7011.47 

20.61-25.65 

27.0711.60 

24.23-30.29 
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